Oh! 



TMUP-HEL-9503 
March 29, 1995 
Revised: June 16, 1995 
hep-ph/ 9504222 



Neutrino Flavor Mixing Constrained by Accelerator and Reactor 
0\ • Experiments and Solar Neutrino Observation* 

a^ 

^ ' Hisakazu Minakata ^ 

1—5 ! 

' Department of Physics, Tokyo Metropolitan University, 

(N ■ 

I 1-1 Minami-Osawa, Hachioji, Tokyo 192-03 Japan 

> 

(N 
(N 

^ , 

^ Abstract 



We derive the constraints imposed on neutrino masses and mixing angles 
by performing a combined analysis of the data from the Los Alamos and 
the other terrestrial neutrino oscillation experiments with the assumption of 
■ the flavor-mixing solutions to the solar neutrino problem. In a three-flavor 

mixing scheme which ignores the possibility of sterile neutrinos, we obtain 
severe constraints on the pattern of masses and flavor mixing of neutrinos. 
For example, we show that in the standard Kobayashi-Maskawa type mixing 
matrix <10~^ independent of the choice of the solar neutrino solutions. 
The constraint from the double 13 decay is also discussed. 
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Recently the Liquid Scintillator Neutrino Detector(LSND) group in Los Alamos has 
reported positive evidence for the neutrino oscillation events z7^ z7g If confirmed 

by continuing runs, it provides the first direct evidence for nonvanishing masses and flavor 
mixings of neutrinos. It is certainly an urgent problem to illuminate all possible consequences 
of this discovery ^ . 

Before hearing this exciting report we have been speculating about the neutrino masses 
and mixing with the help of information of astrophysical neutrinos detected by the various 
underground detectors. Among them the most popular one is the solar neutrino deficit 
which lasts more than 20 years 0], but still refuses resolution by the standard solar model 
0. Now the solar neutrino deficit is confirmed to exist by the four different experiments, the 
chlorine 0, the Kamiokande and the two gallium experiments P,^. A less popular but 
an equally important piece of information comes from the atmospheric neutrino anomaly, a 
deficit in the ratio (i/^ + i7^)/{ve + J^e), which was first claimed to exist by the Kamiokande 
group |]10[, and subsequently supported by the 1MB ||rT| and the Soudan 2 ||12[ experiments. 



The new "multi-GeV" data sample reported by the Kamiokande group provides strong 
support for the interpretation of the atmospheric neutrino anomaly in terms of the neutrino 
oscillation. 

In this paper we discuss possible implications of the LSND result by performing a com- 
bined analysis of the data with that of solar neutrino observation and the reactor exper- 
iments. Our analysis is based on the three-flavor mixing scheme of neutrinos without in- 
troducing any sterile species. We will recognize, as the analysis proceeds, that it is crucial 
to include the information of the disappearance experiments. The most stringent limits 



achieved for z/e and z/^ channels are the ones reported by the Bugey [Tj] and the CDHS |T5 
groups, respectively. 

The most disturbing feature of this restricted framework is, of course, that it cannot 
accomodate all of them, the LSND result, the solar neutrino deficit, and the atmospheric 
neutrino anomaly. The reason is very simple; three mass scales (more precisely, mass-squared 
differences) which are involved in these three phenomena are too far apart from each other 



to be accomodated in the three-generation scheme 0. In the present paper we therefore 
omit atmospheric neutrino data from the analysis. A combined analysis of the alternative 
combination, LSND- Atmospheric neutrino, is presented elsewhere . 

To orient our analysis and to focus some important points we summarize below the 
features of the LSND and the disappearance experiments, and the flavor mixing solution to 
the solar neutrino problem. 
(1) The LSND experiment 



By using low energy beam from stopped muons the LSND group is able to perform 
the appearance experiment — > z7g which can probe the neutrino oscillation parameters up 
to Am^ = 0.2eV^(0.5eV^) for sin^ 26 = 0.01(0.001). Such great sensitivity has been achieved 
by the intense pion beam from the LAMPF accelerator and by the low enough energies of 
z7^ beam. After oscillating into z7g, it produces positrons via the reaction u^p e^n which 
are measured in the energy window 36MeV < Eg < 60 MeV. The produced neutrons are 
captured by protons through the reaction np d7, whose gamma rays of energy 2.2 MeV 
can serve for delayed coincidence. The experimental group may have observed excess in 
gamma ray-correlated events over the estimated background. 
(2) The disappearance experiments [TIJTH]: 



This type of experiment measures the attenuation of neutrino beam from a reactor or 
an accelerator. The intensity of z7g beam from the reactor at Bugey is measured by the 
^Li-loaded liquid scintillator located at 15, 40 and 95m from the reactor [|14|. From the 
viewpoint of neutrino mixing the experiment measures 1 — P{ve i'e)-, where P^C'e ^ i'e) 
implies the survival probability of electron antineutrinos. The resulting bound is rather 
severe and may be summarized as 1 — P^Ue ^ t'e) < 5 % including the statistical and the 
systematic uncertainties. A comparable limit is achieved for channel by the CDHS group 
p!5| using z/^ beam from the CERN PS. 

(3) The flavor-mixing solutions to the solar neutrino problem: 

The most popular and perhaps most appealing solution to the solar neutrino prob- 
lem is the one provided by the Mikhyev-Smirnov-Wolfenstein(MSW) mechanism fl^. It 



utilizes the (by now familiar) mechanism of resonant enhancement of neutrino flavor con- 
version in solar matter. At the moment there exist two options in this type of solution; 
the small-angle (nonadiabatic), and the large-angle solutions. The characteristic values of 
the neutrino mixing parameters are determined by, for example, an extensive analysis by 
Hata and Langacker and are given as Am^ ~ 6 x 10~^eV^, sin^2^^ ~ 7 x 10~^, and 
Am^ ^ 9 X lO^^eV^, sin^ 26 ~ 0.6, for the small-angle and the large-angle solutions, 
respectively. Their analysis is done with the two-flavor mixing scheme, and Am? and 9 
indicate the mass-squared difference and the mixing angle, respectively. We should mention 
that there still exists the possibility of vacuum neutrino oscillation as a mechanism for the 
solar neutrino deficit. This is the old solution 0] but one still alive . 



One of the most important features in our combined analysis of the LSND and the 
solar neutrino solution is that a huge mass hierarchy is involved. (In fact, it is not really 
the hierarchy in mass but is the hierarchy in the mass-difference.) One may classify the 
hierarchy of neutrino masses into two types: 

(a) m\ ~ m\ ^ m\ (b) m\ 3> m\ ~ m\ (1) 

Here the symbols ~ and ^ imply the difference by < lO^^eV^ and ~l-10eV^, respectively. 
In (|I]) we have chosen the third state as the grossly departed mass eigenstate. The other 
types of mass hierarchies which can be obtained by permuting 1,2, and 3 can be taken care 
of by an appropriate choice of angles because they merely represent relabeling of the mass 
eigenstates. Unlike the case of vacuum neutrino oscillation the relative magnitude of 



the masses connected by ~ does have important meaning and will be discussed below. 

We briefly review the neutrino oscillation with three flavors of neutrinos. We define the 
neutrino mixing matrix U which relates the flavor eigenstate //^(a = e,/i,r) and the mass- 
eigenstate = 1,2,3) in vacuum as z/^ = Uaii^i- In this paper we assume CP invariance. 
The evolution equation of the flavor eigenstate takes the form 



where ae(x) = \/2GFNe(x) indicates the effect of matter potential which affects only electron 
neutrinos. Here, Gf,E and A^e denote, in order, the Fermi constant, the neutrino energy, 
and the electron number density in the sun. 

Without matter effect the equation can be easily integrated to yield the oscillation 
probability z/q, — > as 



j>i 



AE 



(3) 



The mass-squared difference Am|,- is defined as Am^^ = m^- — mf{j > i) in this paper. Notice 
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that when we discuss the neutrino evolution in matter the sign of Arn^j does have physical 
meaning, unlike the case of vacuum neutrino oscillation (^. 

With matter effect the equation becomes complicated but it is tractable because of the 
mass hierarchy involved in our analysis. To make our discussion transparent we specify the 
form of the mixing matrix in a form U = U23U13U12, 



(4) 



where U23, f/13, and U12 denote the three matrices in (^), in order, from left to right. Here 
Cij and Sij are the short-hand notations for cos 6 ij and sin 6ij, respectively. This is nothing 
but the standard form of the Kobayashi-Maskawa matrix |2^, which is now employed as 
the neutrino mixing matrix. We note that CP invariance renders the three angles real and 
they can all be made to lie in the first quadrant by an appropriate redefinition of neutrino 
phases. 

The definition of the mixing matrix (^) is convenient in dealing with the mass patterns 
(0) in which 1-2 level crossing is responsible for the solar neutrino deficit. If one wants 
to discuss the other type of mass pattern which can be obtained by permuting 1, 2, and 
3 one may make a different choice of the U matrix (i.e., redefinition of angles) convenient 



for them. For example, U = U13U12U23 for 2-3 level crossing. We note that one can take 
^12 < f so that Ue is close to z/i under the two-level crossing approximation. It does not hurt 
the generality of our analysis because the constraints from the terrestrial experiments to be 
discussed below do not involve Ou- 

We examine the MSW effect in the three-flavor mixing scheme. We multiply, following 
Kuo and Pantaleone |2T|, f/i3^f/23^ the equation @ to obtain the evolution equation for 
the modified neutrino basis = (f^i3^f^23^)/3at'a: 



dx 
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where we have used simplified notation A _ 
the sign of A has physical significance. 



(5) 



and E = |(Am?3 + Am^3). Note that 



From @ one realizes that the effective two-level crossing approximation is justified unless 
sin26'i3ae is extraordinarily large compared with other elements, which is not the case in our 
problem. Moreover, an evaluation of the perturbative corrections to the energy eigenvalues 
due to this off-diagonal term reveals that they are of the order of — (sin 26'i3ae)^/| J2 I 121| 
which is negligible for the mass hierarchy Amf^ ^ Am23 ^ Amfg- Notice that this is 
true for both of the types-a and -b mass patterns given in (|1]). Also it can be shown that 
the correction to the difference between two energy eigenvalues in matter vanishes at the 
resonance point. Therefore, the off-diagonal terms affect the discussion of the adiabaticity 
condition in the effective two-level problem only through higher-order corrections. 

Having established the validity of the effective two-level approximation we proceed to 
the combined analysis. We first derive the approximate formulas for the terrestrial neutrino 
experiments, taking into account the mass hierarchy and the experimental parameters. With 



mass hierarchy (0) the oscillation probability consists of two terms, one involving large A 



m 



13 



and other small Amf2- The coefficient of the former term have a simplified expression due 
to the orthogonality of the mixing matrix [|T^. The latter term is smaller by factors of 
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(Am^2/^^i3)^ — 10~^° compared with the former. So it can be safely ignored. 

In certain cases, the argument of the sine function with larger mass difference takes the 
large values, e. g., 10-100 for Am^g = 1 — lOeV^ with E= 4MeV and L= 40m, the typical 
parameters in the Bugey experiment. Therefore, it can be averaged to be | and we obtain, 
as the formula for the Bugey disappearance experiment, 

1 - P{v, -> Z7,) = 2C?3S?3. (6) 

For the LSND experiments the arguments of sine factors with Amfg and Am^a are of order 
unity and so we cannot average. The oscillation probability in the LSND experiment is thus 

Similarly, the oscillation probabilities of and Vr channels take the forms 

1 - P{y, - y,) = 444(1 - 44) sin2(^^), (8) 

P{v, V,) = 4444 sin2(^^). (9) 
From the data of the Bugey and the LSND experiments described earlier we require 

44 < 2.5 X 10"^ = 5 (10) 

444 = e<io-^ (11) 

where we have defined 5 and e such that 25 and 4e imply the attenuation of z7g in the Bugey 
and the "rate" in the LSND experiments, respectively. Since the latter quantity is still 
subject to the uncertainty p8| , p9| , we take a conservative attitude and assume that e is less 
than or equal to ~ 10~^. From (|1^) and (|11]) it follows that ^13 must be either small or close 



tof. 
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e < 



■^13 
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< 2.6 X 10^2 ~ 5 (12) 



The angle S23 is severely constrained by the the disappearance experiment. The 
CDHS group [|T^ obtained the constraint on sin^ 2^ (in two-flavor scheme), sin^20<O.l in 
the mass range leV"^ < Am^g < lOeV^, where the constraint from the CDHS experiment is 
most stringent. Using this with (|^) we obtain the bound on S23 as 

\ 

(13) 
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In large- Am^3 region a better bound can be obtained from z/„ — >■ z/^ oscillation experiment. 



Using the data of Fermilab E531 experiment with (^ we obtain an approximate expres- 
sion of the Am^-dependent bound in the same mass range: 



2 „2 „4 ^noc / ^"^13 
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(14) 



We will show below that the solar neutrino solutions force us to choose the small-Si3 



solution out of (|T2|). Using this information apriori with (O) and (|I4|) we obtain the bound 



on s 



23' 



f.2 \ 



'23 



<2.5 X 10" 



y ^23 J 



\3.i 



3eV^ 



(leV^ < Amfg < S.SeV^) 
(3.3eV^ < Am^g < lOeV^ 



(15) 



Now let us turn to the constraint implied by the solar neutrino solutions. Using the local 
two-level crossing approximation established above Am^2 ^"^^ ^12 determined to be 



Amf2 ^ 6 X IQ-'^eV^ s 



12 



1.75 X 10" 



(small- angle) 



(16) 



and 



^12 



9 X IQ-'^eV^ ^ 



0.184, 



(large- angle) 



(17) 



for the small-angle and the large-angle MSW solutions, respectively. Of course, the mass- 
squared difference Am^2 must be positive in order that the MSW mechanism acts for neu- 
trinos, not for antineutrinos. This implies > ml between almost degenerate neutrino 
states in (|^). 



So far we have discussed the constraints from the terrestrial experiments and that from 
the MSW solar neutrino solutions separately (except what we did for S23). While we have 
obtained highly nontrivial constraints (|1^)-(|T^), it would be more interesting if there arise 
further restrictions by considering these three experimental requirements simultaneously. 
We show that this in fact occurs. 

We demonstrate that the large-si3 solution (i. e., cf^ < 2.6 x 10~^) is not acceptable as 
a solution to the solar neutrino problem. To show this we note the relationship between the 
flavor basis Ua and the basis i7a introduced as a convenient basis for describing the effectively 
local two-level resonance in neutrino evolution in matter. It is 
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(18) 



It we take the Iarge-Si3 solution the initial condition [Uf. v^-] [1,0,0] at the solar 

core is translated into the initial condition [ue z7^ i^t] = [v^, 0,1]. Since effectively 
decouples with 2x2 submatrix which has resonance, as can be seen in (^, it experiences 
no significant change in the solar interior. Therefore, z/g's which departed the solar core 
just leave the sun with no appreciable attenuation in their fiux.Q So the Iarge-Si3 solution 
cannot explain the solar neutrino deficit. This conclusion is valid for both of the small-angle 
and the large-angle MSW solutions, because the relation (|l^) is independent of S12. Thus, 
we are left with the small-si3 solution with the additional constraints ([TsD together with 
either (|1^), or (|1^, depending upon the small-angle and the large-angle MSW solutions, 
respectively. 



^One can reach the same conclusion by using the local two- level crossing representation of 1-3 
channel |21|. In this case "^v" does not experience a level crossing not because it decouples but 
because the resonance condition cannot be met due to the mass hierarchy Am^2 ^ ^"^13 ^ ^^23 
with which we are working. 
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Now we turn to the discussion of the solar neutrino solution based on the vacuum neutrino 
oscillations. Because of the mass hierarchy involved, 1-lOeV^ ~ ^"^13 ~ ^"^23 ^ ^^12 — 
lO^^'^eV^, the relevant oscillation probability can be written as 

P{^e - ^e) = 1 - 2cl,sl, - 4cL^?24sin^(^^), (19) 

where the sine-squared terms with Am^g and Am^g are averaged. It is an excel- 
lent approximation because the argument of sine takes a very large value, ~ 2 x 
10^°(AmVleV^)(L/lAU)(£'/10MeV)~\ In view of ([ID and noticing the Bugey constraint 
(|T0|) we realize that the only possible way of having solar neutrino deficit of ~50% level is 
to have an effective two-fiavor description of P{i'e '^e)- Namely, we have to demand 

c^2Si2 — i(= maximum value), 

•-IS — -L) 



to have a gross deficit. Under the condition (|20|) the vacuum neutrino oscillation (O) will 



provide an acceptable solution to the solar neutrino problem (with possible slight changes 



in mixing parameters) as shown in [^. Thus, the combined analysis of the LSND and the 
Bugey experiments with the vacuum mixing solution of the solar neutrino problem again 
prefers the small-sis solution. 

An additional constraint arises in the case of Majorana neutrinos. The quantity 

3 

< m^e > =Y1 VjlUejl'^rrij 
i=i 

= '7iC?2C?3"^i + ^2S?2C?3m2 + r]^s\^m^ (21) 

is constrained to be less than ~1 eV by the non-observation of the neutrinoless double (5 
decay in various experiments Notice that we are working with the representation in 



which the mixing matrix is real under the assumption of CP invariance, and r^j = ±1 in (^11) 
is the CP phase. 

The constraint from the double (3 decay acts differently for the type-a and the type-b 
mass hierarchies in (^. In the type-a case there is a chance for cancellation between nearly 

10 



degenerate two masses, but there is no chance in the type-b case because the heavy mass is 
carried by a unique mass eigenstate. 

We, however, encounter a new situation in the consistent solutions obtained in our com- 
bined analysis. For the type-b mass hierarchy the double (3 decay constraint is automatically 
satisfied by the small-s^g. On the other hand, the constraint for the type-a mass hierarchy 
has nontrivial consequences. In the case of small-angle MSW solution, the cancellation be- 
tween 1 and 2 mass eigenstates is hopeless because ^'^'^ tiny; — 10"^. In the case of 
large-angle MSW solution the situation is better but we still obtain < uij^f. >~ 1.7 eV for 



mi = m2 = 2.4 eV, which is larger by factor of ~ 2 than the experimental bound p4[. In 
the case of vacuum oscillation solution, there exists better chance for cancellation because 
is large. The double (3 constraint can be cleared with a mild condition 33° < 612 < 57°. 
Thus, the double j3 constraint prefers the vacuum oscillation solution for the type-a mass 
hierarchy, while it is automatically satisfied for the type-b mass pattern. 

In this paper we have discussed the constraints imposed on neutrino masses and mixings 
when we demand the consistency with the LSND and the other terrestrial neutrino exper- 
iments and the fiavor mixing solutions to the solar neutrino problem. Independent of the 
choice of the three solar neutrino solutions, the small- and the large-angle MSW, and the 
vacuum oscillation solutions, is constrained to be small, s\^<5 ^ 10^^, for different rea- 

""""^^ ''23 



sons in the MSW and in the vacuum solutions, respectively. The angle is subject to the 



constraint (^) which says that it is either small, < 10^ or large, ~ 1. On the other hand, 
the value of depends upon the solar neutrino solutions; — 1-8 x 10^^, 0.18, 0.5 for the 
small-angle MSW, the large-angle MSW, and the vacuum oscillation solutions, respectively. 

The physical interpretation of the solution is clearest in the small-angle MSW solution. 
The vi state is approximately identical with "z/g" , and z/2 and z/3 are in general mixtures of 
'V^"and "z/^". Therefore, the type-a contains an inverted mass hierarchy and the type-b 
implies a normal mass pattern. Likewise one can obtain analogous physical interpretations 
in other solutions but the vi states is less pure with larger mixing angles. 

Finally, a few remarks are in order: 
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(1) The present analysis is less powerful in constraining the absolute values of masses than 
the relative masses and the mixing angles, a general feature noticed in . The constraints 
obtained in this paper do allow, for Dirac neutrinos, the type-b solution with e.g., mi = 6 
eV, m2 = 6 eV, and ma = 6.5 eV, which is consistent with the direct mass measurement 

and the cosmological constraints [^]. This extreme choice would provide the possibility 
that the light neutrinos could fill the entire part of the missing mass in the universe, but 
with possible troubles with galaxy formation ||26|,|3[|. 

(2) One must be careful in comparing the resulting constraints obtained in this paper with 



that of Ref. not only because they stand on entirely different basis but also because the 
definition of the angles are different. This complication arises due to the fact that we are 



working with the different definitions of the mass eigenstates here and in Ref. [0. We are 



planning to present a unified and more transparent description in the future [27 



I thank Hiroshi Nunokawa for discussions, Eligio Lisi and the referee of Ref. |]T6[ for 
calling my attention to various terrestrial experiments which have not addressed in the 
earlier version of this paper. This work is supported in part by Grant-in-Aid for Sciectific 
Research of the Ministry of Education, Science and Culture, #0560355. 

Note Added: After submitting the earlier version of this paper there appeared two reports 
from the LSND group [^,^] with mutually conflicting conclusions. 
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